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Abstract 

A general nonautonomous Mackey-Glass equation for the regula¬ 
tion of the hematopoiesis with several non-constant delays is studied. 
Using topological degree methods, we prove the existence and multi¬ 
plicity of positive periodic solutions. 


1 Introduction 


The following nonlinear antonomous delay differential equation was proposed 
by Mackey and Glass to study the regulation of hematopoiesis: 


dPjt) xe-Pit - r) 
dt 0- + P-(t-r) ^ 


( 1 ) 


Here X,6,n,'y,T are positive constants, P{t) is the concentration of cells in 
the circulating blood and the flux function f{v) = of cells into the 

blood stream depends on the cell concentration at an earlier time. The delay 
r describes the time between the start of cellular production in the bone 
marrow and the release of mature cells into the blood. It is assumed that 
the cells are lost at a rate proportional to their concentration, namely 'jP{t), 
where 7 is the decay rate. This equation is a model of ‘dynamic disease’. 

*This work was fully supported by the projects UBACyT 20020120100029BA and PIP 
11220090100637 CONICET 


1 






This type of population dynamics equations attracted the interest of many 
researchers. Different aspects and properties of Q have been studied by 
various authors, see for example [5|[^[l0]. 

Most often, the environment is not temporally constant; thus, it is intu¬ 
itive to assume that this fact influences many biological dynamical systems 
and suggests the need of considering time-dependent parameters. Moreover, 
as remarked in IdllS, 16 , more realistic models are those in which periodicity 


of the environment and time delay play a role (for more details, see e.g. [l^). 
In view of this, the following model is proposed in [TB]: 


x'{t) = 


q{t)x{t) 


r + x'^{t — mT) 


-p{t)x{t) 


( 2 ) 


where m, n are positive integers, p, q are positive T—periodic functions and 
the delay r := mT is a multiple of the period determined by the environment. 

In order to establish a more realistic model, it is convenient to introduce 
a more general delay that extends the two above-referred cases. Instead of 
assuming that the delay is constant or a multiple of the period of the envi¬ 
ronment, more general models are obtained by assuming that the time delay 
r is an arbitrary continuous nonnegative T—periodic function depending on 
t. The more general equation 


x'{t) 


a{t)x{t — T{t)) 
1 - 1 - x'^(t — r(f)) 


b{t)x{t) 


( 3 ) 


where a, b and r are continuous positive T—periodic functions was studied for 
example in 17 21 . Different aspects of equation ([^ have been considered; 


in particular, existence of positive T-periodic solutions was proven, in most 
cases, using appropriate fixed point theorems. In 20 , coincidence degree 


theory was employed to prove the existence of a positive T—periodic solution 
under a condition that can be regarded as a particular application of Theorem 
2.2, case (2) below. Moreover, when a{t) = 'yb{t) for some 7 > 0 and when 


r, a and b are constant, the conditions 7 > 1 and a > b respectively are both 
necessary and sufficient for the existence of positive T-periodic solutions. 
The following more general model was studied in and [^: 


M 


x'{t) = 


k=l 


rk{t)x^{t - gk{t)) 
1 + x'rit - gk{t)) 


b{t)x{t). 


( 4 ) 
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Here, 7 is a positive constant and rk,b are positive T—periodic continuous 
functions. For 5 = 1, existence and uniqueness of positive T-periodic so¬ 
lutions was studied in for the particular case of constant proportional 
delays Qk = hT] moreover, for general continuous, positive T—periodic gk, 
attractiveness of some specific positive periodic solutions was studied. For 
the case 5 = 0 and continuous positive and T—periodic, existence and 
uniqueness of positive T—periodic solutions of Q was proven in by fixed 
point methods, provided that one of the following conditions is satisfied: 


b(u)du pT ^ \ ~ 

.T° 

e/o b{u)du _ ^ ^ J 

Motivated by the preceding discussion, we consider the following more 
general nonlinear nonautonomous model with several delays 

// \ /\ (t — Tkit)) ,/X / N 

^ = E hn(t)-b(t)Mt) (5) 

where rk(t),b{t),Tk(t) and /ifc(t) are positive and T-periodic functions and 
Xk,mk,nk are positive constants. 

Existence of solutions of (|^ under appropriate conditions follows from 
several abstract results, although multiplicity results are more scarce. For 
example, in and Krasnoselskii type hxed point theorem in cones 
were employed in order to obtain conditions for the existence of at least 
two T—periodic solutions of the general equation 

x'{t) = -a{t)x{t) + f{t, x{t - ...,x{t- (6) 


It is observed, however, that these results can be applied only to few par¬ 
ticular sub-cases of ([^. In such cases, the conclusions are comparable to 
our results below. Moreover, the existence of three nonnegative periodic 
solutions of ([^ was studied by using Leggett-Williams fixed point theorem 
in 


15 . However, the conditions obtained in 13 , as pointed out by the 


authors, are very difficult to apply to ^ with M = l,m = l,r = p. Thus, 
they established a complementary result with more straightforward condi¬ 
tions that can be applied to this model. Unfortunately, in 14 
observed that this latter result was incorrect. 


the authors 
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In section 14, Applications], the hematopoiesis model ([^ for M = 1, r = 
/i was studied, and the conditions obtained by the authors are similar to 
the ones proposed in Theorem 4.2 (1) below. It is worth mentioning that. 


in the referred result of 14 , only two of the three T-periodic solutions are 


positive and the third one is positive if /(t,0) is not identically zero. This 
assumption is very restrictive and clearly is not fulhlled in ([^. Moreover, 
all the mentioned works do not contemplate the superlinear case of ([^ (that 
is, irik > Uk + 1 for some k). Finally, we may mention the results in |^ , 
in which the existence of at least 2n solutions of ([^ is proven, although the 
conditions are not applicable to our model. 

Our goal in this paper is to establish sufficient criteria to guarantee, on 
the one hand, the existence of positive T-periodic solutions of ([^ and, on 
the other hand, multiplicity of such solutions. Using degree theory, we shall 
obtain a set of natural and easy-to-verify conditions for the existence of one 
or more solutions. 

The following notation will be used throughout the paper. Let 

Ct '■= {uit) G ^(M, M) : u{t + T) = u{t) for all t} 

denote the space of continuous T—periodic functions and define, for r < s, 

:= {u{t) E Ct ■ r < u(t) < s for all t} . 

The closure of shall be denoted by d{X^). The average, the maximum 
value and the minimum value of an arbitrary function ip E Ct shall be 
denoted respectively by Pmax and Pmim namely 


(p := — p{t) dt, Prnax = max(y9(t), Prnin = min9?(t). 


T 


/o [OT] [0,r] 


In order to simplify some computations, we set y{t) := ln(a;(t)) and trans¬ 
form ([^ into the equivalent equation 


M 


y{t) = 'Y^^krk{t) 


k=l 


X -)- pri^vd-k'kd)) 




(7) 


Finally we define, for convenience, the function 0 : M —>■ M by 

M 


0(7) := ^ y^kTk— —V- - h. 

k=l 


' 1 -f 


( 8 ) 
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The proof of our results shall be based on the continuation method. 
Specihcally, we shall apply the following existence theorem, established in . 

Theorem 1.1 Assume there exist constants 71 < 72 such that 

1. If y E cl{Xff) satisfies 

f grnky{t-Tkit)}-y(t) \ 

»'(*) “ 1 + en.l,(.-»C)) “ ''(*> j 

for some a G (0,1], then y G Xff. 

2 . 0(71)0(72) < 0 . 

Then has at least one solution in Xff. 

Roughly speaking, if 0 has different signs at both ends of some interval 
[ 71 , 72 ] C M then the continuation theorem guarantees the existence of a T- 
periodic solution 1 / of ([^ such that y{t) G ( 71 , 72 ) for all t. However, the hrst 
condition of Theorem |1.1| requires, in some sense, that the sign of 0 does not 
change too fast. 

The main part of our analysis shall be based on a study of the behavior 
of 0. For the existence of solutions it suffices, in most cases, to consider its 
behavior at ±cxd; for the multiplicity results a more careful study is needed, 
in order to hnd “large enough” intervals of positivity and negativity of 0. 
With that end in mind, we shall consider the sets 


Ml := {A; : 0 < mfc < 1}, M 2 := {k : mk = 1}, M 3 := {k : 1 < mk < Hk + 1} 


M4 := {k : mk = nk + 1}, M5 := {k : mk > Uk + 1} 


and 

g(mfe-l)7 

^.(7) := 7 ^ 

k&Mi ^ 

SO we may write 0(7) = For notation convenience, we also 

dehne C := Tb = b(t) dt. 

This setting proves to be useful, since the limits lim..y_,.±oo 007) are easy 
to compute and, moreover, 007 ) is strictly monotone for *03 and a sum of 
one-hump functions for i = 3. Thus, the behavior of 0 can be understood by 
studying the interaction of these different terms. 
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The paper is organized as follows. In the next section, we adapt the 
abstract continnation theorem |[^ Thm 2.1] to equation ([^ and prove the 
existence of positive T-periodic solutions for the different cases. In Section 
we give sufficient conditions for the existence of 2, 3 or 4 positive T-periodic 
solutions. Finally, in Section]^ we give an example with at least 6 positive 
T-periodic solutions. 

2 Existence of positive T—periodic solutions. 

In order to present our existence results in a more comprehensive way, we 
shall consider three different cases: the superlinear case (m^ > Uk + l for some 
k ), the the sublinear case {rrik < rik + 1 for all k) and the asymptotically 
linear case (mfc < rik + 1 for all k and mj = rij + 1 for some j). We give a 
detailed proof only of the hrst result, since the other two follow similarly. 

Theorem 2.1 Assume mj > rij + 1 for some j. Furthermore, assume 
that one of the following conditions is fulfilled: 

1. mk > 1 for all k. 

2. mk > 1 for all k, mi = 1 for some i and J 2 keM 2 < b{t) for 

all t. 

3. mi < 1 for some i and Yl!k=i < bit) for all t and 

some constant 71 . 

Then 0 admits at least one positive T—periodic solution. 

Theorem 2.2 Assume mk < Uk + l for all k. Furthermore, assume that 
one of the following conditions is fulfilled: 

1. mi <1 for some i. 

2 . mk > 1 for all k, mi = 1 for some i and J 2 k&M 2 ^krkft) > bft) for all 

t. 

3. mk > 1 for all k and 

g(mfe-l)7i 

for all t and some arbitrary constant 71 . 
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Then admits at least one positive T—periodic solution. 

Theorem 2.3 Assume + 1 for all k and mj = nj + 1 for some 

j. Furthermore, assume that one of the following conditions is fulfilled: 


1. mk > 1 for all k and ^krk{t)e > b{t) for all t. 

2. mk > 1 for all k, mi = 1 for some i, Yhk&M^ ^kf'k{t)e~^'^'^ > h{t) and 
EfceMa >^krk{t)e^ < b{t) for all t. 

3. 0 < mj < 1 for some i and ^k^k(t)e‘"^'^ < b(t) for all t. 


Then & admits at least one positive T—periodic solution. 

Proof of Theorem 2.1\ Let ?/ be a T—periodic solution of (|^ with 0 < cr < 1, 
then y'{t) > —b{t) and hence y{ti)—y{t 2 ) < b{t)dt for any ti <t 2 < F+T. 

This implies, since y{t) is T—periodic, that ymax — l/min < Jo Ht)dt = C. 
Moreover, since mk > Uk + 1 for some k it follows that 0 ( 7 ) > 0 when 7 is 
large enough. Assume that ymax is achieved at some value t*, then 




2_^^k^kit ) 


and consequently 


> 


M 

y^^hrkjt*) 

k=l 


^rrikiy max c) 
-j— ^Tlj^ymax 


p{mk-l)ymax p-Cm^ 

bin>j2^krk{n —. 

r J _L pitjf^ymax 

k=l 

Again, since m^ > n^ + l for some k we deduce that ymax cannot be too large. 
Thus, we may £x 72 3> 0 such that ymax < 72 for every y G Ct satisfying @ 
and 0 ( 72 ) >0. In a similar fashion, we look for 71 < 72 such that 0 ( 71 ) < 0 
and ymin 7 ^ 7 i- 


Case 1: mk > 1 for all k. Here 

0 ( 7 ) —)■ —6 as 7 —)■ — 00 . 
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Let y G Ct be a solution of ([^ and fix f* such that y{t^) = ymin, then 


M 


bit*) < X] hr kit* 


k=l 


Suppose that ymin = 7i, then 


M 


bit*) < Afcrfc(t) 


k=l 


—|— ^^kymin 


g(mfe-l)7i+Cmfc 

1 + 


The right-hand side of the latter inequality tends to zero as 71 —)■ 
—00. We deduce that ymin cannot take arbitrarily large negative values; 
hence, it suffices to take 71 -C 0 . 

Case 2 . > 1 for all k and ruj = 1 for some j. In this case, 

0(7) -t hh -b < hhe^ -b <0 

fceM2 keM2 

7 —>■ —00. On the other hand, if y G Ct satishes 0 then 


as 


bit*) > Y. = yj hr,iu 


nC 


fceM2 


keM2 


X -\- gnfe(ymin + C) 


and, again, we deduce that ymin cannot take too large negative values. 
Thus, it suffices to take 71 -C 0 . 

Case 3 . rok < 1 for some k. From the hypothesis. 


M 


g(mfe-l)7i _ g(mfc-l)7igmfcC 

0 ( 71 ) = 2 ^ - 7 — -b<Y hr-k— -& < 0 . 


■ 1 -I- gnfc7i 
k=l k=l 


1 - 1 - 

Moreover, if ymin is achieved at some value t*, then 


M 


_^ g(^fc ^)ymin 

bit*) < —Y 


k=l 




We conclude that ymin 7^ 7 i- 









rem 

2 ’. 


Remark 2.1 It is easy to verify that that the second condition in Theo- 

< b{t) and 


2.3 can he replaced by 


rrik > 1 for all k,mi = 1 for some i, Y.k&M 4 
EfceMa ^kTkit) > bit) for all t. 


3 Multiplicity 


In this section, we shall employ Theorem 1.1 to prove the existence of multiple 


solutions. It is worth noticing that, when 0 is monotone, it changes sign at 
most once and the method cannot be applied. When 0 is not monotone, it 
is not enough to obtain intervals of positivity and negativity: as mentioned, 
it is required that 0 does not change sign too rapidly. For a more detailed 
analysis, the following functions shall be helpful: 


«(7,^) ■.= 2^Xkrk{t) 

JtL pimk-^hpCrrik 

/3h,t) :=^Afcrfc(t) -b{t) 

As before, our results shall be split in three different theorems, for the 
superlinear, sublinear and asymptotically linear cases. More concretely. 


Theorem 3.1 Assume that mj > Uj + 1 for some j. 

1. Let mk > 1 for all k and 1 < m* < Uj + 1 for some i. Assume there 
exist constants 71 < 72 such that 

0 ( 71 , f) > 0 > 0 ( 72 , t) for all t. 

Then admits at least 3 positive T—periodic solutions. 

2. Let mk > 1 for all k, = 1 for some i, mk 0 (1,'n.fc + 1) for all k. 
Assume that 

XkTkif) > h{t) for all t 

fceM2 
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and there exists 71 such that 

/ 3 ( 7 i, t) < 0 for all t. 

Then 10 admits at least 2 positive T—periodic solutions. 

3. Let ruk > 1 for all k, rui = 1 for some i and 1 < <71^ + 1 for some 

s. Assume that 

Xkrk{t)e^ < b{t) for all t 

k&M2 

and there exist constants 71 < 72 such that 

Q!( 7 i,t) > 0 > / 3 ( 72 ,t) for all t. 

Then 0 admits at least 3 positive T—periodic solutions. 

4 . Let mi < 1 for some i, mk ^ (1,^^ + 1) for all k. Assume 

/ 3 ( 7 i,t) < 0 for all t and some constant 71 . 

Then admits at least 2 positive T—periodic solutions. 

5. Let mi < 1 for some i, 1 < < 77 ^ + 1 for some s. Assume there exist 

some constants 71 < 72 < 73 such that 

(^{ 32 , t) > 0 > /9(7j,t) for all t, 

for i = 1,3. Then admits at least 4 positive T—periodic solutions. 

Theorem 3.2 Assume that mk <77^ + 1 for all k. 

1. Let mk > 1 for all k and assume there exists a constant 71 such that 

a( 7 i, t) > 0 for all t. 

Then admits at least 2 positive T—periodic solutions. 
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2. Let ruk > 1 for all k, rrii = 1, nij > 1 for some i,j. Assume that 

< h{t) for all t 

keM2 

and there exists a constant 71 such that 

a( 7 i, t) > 0 for all t. 

Then admits at least 2 positive T—periodic solutions. 

3. Let 0 < mi < 1, mj > 1 for some i,j. Assume there exist some 
constants 71 < 72 such that 

(^{ 32 , t) > 0 > / 3 ( 7 i,t) for all t. 

Then admits at least 3 solutions. 

Theorem 3.3 Assume that < Uk + 1 for all k and mj = nj + 1 for 
some j. 

1. Let mk > 1 for all k and 1 < m* < tt.* + 1 for some i. Assume that 

< b{t) for all t 

keMi 

and there exists a constant 71 such that 

0 ( 71 , t) > 0 for all t. 

Then admits at least 2 positive T—periodic solutions. 

2. Let 0 < mi < 1 for some i, mk ^ (1, Uk + 1) for all k. Assume that 

> hit) for all t 

keM4 

and there exists 71 such that 

/ 3 ( 7 i, t) < 0 for all t. 

Then 0 admits at least 2 positive T—periodic solutions. 
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3. Let 0 < m* < 1 and 1 < <71^ + 1 for some i, s. Assume that 

< b{t) 

fceM4 

and there exist constants 71 < 72 such that 

tt( 72 ,^) > 0 > / 9 ( 7 i,t) for all t. 


Then has at least 3 positive T—periodic solutions. 


As before, we shall only prove the hrst case of Theorem 3T, since all the 
remaining cases follow in an analogous way. 

Proof of Theorem 5A, case 1\ We shall apply Theorem 11.11 on open bounded 


sets X3f and Xf^, with 70 < 71 and 73 > 72 to be determined, 
begin, observe that 

0 ( 7 ) —)■ — 6 as 7 —)■ —00 

and 


To 


0 ( 7 ) —>■ +CXD as 7 —)■ + 00 . 


In the same way of Theorem 2.1 it is proven that, if 70 ^ 0 then there 
exists y G Xf^^ solution of ([7). 

On the other hand, for ml t it is seen that 


0(71) > a( 7 i,t) > 0 . 

Moreover, ii y E cl{Xf^) is a solution of ([^ with 0 < a < 1 and ymin = y{t*)-i 
then 


M 


= '^XkrkiU 


M 


> 


^ ^ Xk'l'kif* 


k=l 


o'^kymin 


^m^y(U-Tk{U)) 


M 


k=l 


X Q-n-kiy min +c) 


> 


y^^Xkrkjp) 


k=l 


^PPkymin ^ Cuif^ 

X ^rikiy min +c)- 


It follows that ymin 7 ^ 7 i- 
Furthermore, 


0 ( 72 ) < /5(72,^) < 0 
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for all t and we deduce as before that Umax 7 ^ 72 - 

Finally, the existence of 73 S> 0 such that the problem has a solution 
y G follows as in Theorem 


2.1 


The following lemma shows, in the context of Theorem 3.1 (case 1), that 
if Tfc, rrik and are given, then it is possible to hnd parameters Xk such that 
assumptions are fulhlled. Analogous arguments are valid for the remaining 


cases. 


Lemma 3.1 Letvk, 6 : M —)■ M>o be continuous andT—periodic functions 
and mk,nk G M>o such that ruk > 1 for all k, 1 < ruj < nj + 1 for some j, 
rui > Hi + 1 for some i. Then there exist Xk and 71 < 72 such that 

a( 7 i, t) > 0 > /d( 72 , t) for all t. 


Proof Using the sets Mi as before, we may write a and (5 as 

5 5 

«(7, t) = Y^ ai(7, t) - b{t), /3(7, t) = ^( 7 , t) - b{t). 

i=l i=l 

Observe that, for each t G [0, T] and i = 1,. .., 5, ai{-, t) and /?*(•, t) have the 
same qualitative behavior as 0 *. 

We begin by setting the parameters Xk G M 3 . For arbitrary 71 , take 
Xk G M 3 large enough such that 

g(mfe-l)7ig-C'"Xfc 

«3(7i,^) = 2 ^ Krk{t) - bit) > 0 . 

keMs 

For e G (0, bmin), there exists R> 71 such that 

g(r7n.-l)7 Cmfe 

A(7.«) = 1 + e".h-c) ^ 1 + e-.(T-c) < ‘ 

fceMa fceMs 

for 7 > R and all t. Thus, we may £x 72 > i? and proceed with the remaining 
parameters. 

Next, for k G M 4 U M 5 we set Xk small enough so that 

g(m.fe-l)72gCmfe 

2_^ ^kirk)max ^ ^nkp2-C) ^ ~ 

keM4UM5 
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and hence 


g(mfc-l) 72 gCmfe 

(/^4 + /55)(72,t) = 2^ Afcrfc(t) ^ ^ < bmin - 26 < h{t) - 2e. 

keMiUMs 

Tims the conclusion follows since 

/^(72) t) = (/^s + /?4 + /^s)(72, t) — b{t) < e — 2e < 0 

and 

«(7i,f) = («3 + ^4 + a5)(7n^) - K^) 


> 03(71, t) - b{t) > 0 


4 Example 


As shown in Theorem 3.1, case 5, equation (|^ has at least 4 positive T- 


periodic solutions. The following example shows that, in fact, the problem 
may have more solutions. Let k = 4 and b{t) = 1.1 + 0.02 cos(^), T = 0.005, 
mi = 0.95, ni = 2, Airi(t) = 0.04 + 0.002cos(^), m 2 = 4.73, n2 = 3.74, 

ns = 10.2, X3rs{t) = 0.9 + 


X 2 r 2 it) = 1.3 + 0.002cos(^), m3 = 1.0001, 


0.002cos(^), m4 = 1.12, 114 = 0.11, A4r4(t) = 0.06 + 0.002cos(^). 

Set 7i = —5, 72 = —0.3 73 = 0.2, 74 = 5, 75 = 34. It is verihed ( 
Figure 1) that 

0^(72, > 0.09, 0(74, f) > 0.1 for all t 


see 


and 


/9(7i,f) < —0.08,/3(73, f) < —0.01,/3(75, f) < —0.01 for all t. 

Moreover, since 0 < mi = 0.95 < 1 and m4 = 1.12 > 77,4 + ! = 1.11, it follows 
that 

lim 0(7) = lim 0(7) = +00. 

7^—00 7^+00 

Thus, we conclude that (|^ has at least six positive 0.005—periodic solutions 
for arbitrary nonnegative 0.005-periodic delays rfc,/ifc. 
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Figure 1: and I3{'y,t) for each f G [0 : 0.01 : T], 
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